Inclusive cross section and double helicity asymmetry 
for tt" production in p + p collisions at ^/s = 62.4 GeV 
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The PHENIX experiment presents results from the RHIC 2006 run with polarized p-\-p collisions 
at ~ 62.4 GeV, for inclusive tt" production at mid-rapidity. Unpolarized cross section results 
are measured for transverse momenta pr ~ 0.5 to 7 GeV/c. Next-to-leading order perturbative 
quantum chromodynamics calculations are compared with the data, and while the calculations are 
consistent with the measurements, next-to- leading logarithmic corrections improve the agreement. 
Double helicity asymmetries All are presented for pr = 1 to 4 GeV/c and probe the higher range 
of Bjorken x of the gluon (3:3) with better statistical precision than our previous measurements 
at = 200 GeV. These measurements are sensitive to the gluon polarization in the proton for 
0.06 <Xg < 0.4. 

PACS numbers: 13.85.Ni,13.88.+c,21.10.Hw,25.40.Ep 



I. INTRODUCTION 

Spin is a property of particles as fundamental as charge 
and mass. The spin of the proton was first determined 
in the 1920s, yet we still do not have a detailed under- 
standing of what inside the proton makes up the spin of 
the proton. Polarized lepton-nucleon deep-inelastic scat- 
tering (DIS) experiments have revealed that only ^^25% 
of the proton spin can be attributed to the spins of the 
quarks and anti-quarks [l|, Q indicating that the proton 
spin must be largely carried by the spin of the gluons 
and/or orbital angular momentum of quarks and glu- 
ons. Polarized proton-proton collisions at the Relativistic 
Heavy-Ion Collider (RHIC) provide a laboratory to study 
the gluon-spin contribution to the proton spin structure, 
AG, with strongly interacting probes via measurements 
of double helicity asymmetries {All) Hi- 

The All of tt^'s is defined as 



4^ — 



(1) 



where a++{a^ ) represents the 7r° production cross sec- 
tion in polarized p -I- p collisions with the same (opposite) 
helicities. In leading order (LO) perturbative Quantum 
Chromodynamics (pQCD), tt*^ production is the sum of 
all possible subprocesses ab — > cX, where a,b represent 
the initial partons in the protons, c is the final state par- 
ton which fragments into a tt^, and X is the unobserved 
parton. Then All is calculated as 



A'^ — 
^LL — 



-'a.b,c 



ab^cX]~[o.b^cX] 



Df 



(2) 



where fa,b represent unpolarized parton distribution 
functions (PDF) of parton a, b and Afa,b represent po- 
larized PDFs, D^" is a fragmentation function (FF) of 
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parton c to 7r°, o-!'^''^^^] and a^^^'^'^' denote respectively 
the cross section and All of the subprocess ab cX. 
The sum is performed for all possible partons (quarks 
and gluons). The Bjorken- a; dependence of the PDFs, 
the kinematical dependence of the FFs, and the integral 
over all possible kinematics are omitted in the equation. 
The partonic quantities a and ull can be calculated in 
pQCD. Since tt*^ production is dominated by gluon-gluon 
and quark-gluon scattering in the measured pt range 
{pt < 4 GeV/c), All is directly sensitive to the po- 
larized gluon distribution function in the proton. 

Cross section measurements at RHIC have established 
the validity of using a next-to-leading order (NLO) 
pQCD description at y/s — 200 GeV for inclusive mid- 
rapidity 7r° [j, Q and forward tt^ production and 
for mid-rapidity jet Q and direct photon production [1]. 
However, at lower center of mass energy, NLO pQCD 
calculations have been less successful in describing the 
data 0. The inclusion of "threshold resummation" 
at next-to-leading logarithmic accuracy (NLL) [lo| im- 
proves the agreement between theory and data at fixed- 
target energies. While taking into account threshold log- 
arithms at the fixed-target kinematic region is essential, 
they may also need to be accounted for at y/s = 62.4 
GeV, but will provide a smaller effect pl| . 

A precise measurement of the inclusive tt" production 
cross section at ^/s = 62.4 GeV is important for the 
heavy- ion program at RHIC. A new state of dense mat- 
ter is formed in Au-|-Au collisions at 200 GeV and par- 
ton energy loss in the produced dense medium results 
in high px leading hadron suppression. Measurements 
of high Pt data at lower energies are of great impor- 
tance in identifying the energy range at which the sup- 
pression sets in. However, they require solid measure- 
ments of the cross section in p -\- p collisions as a base- 
line for medium effects. At the ISR, inclusive neutral 
and charged pion cross sections were measured several 
times at ~ 62 GeV [l^, [l^ , but they have large un- 
certainties and have a large variation [ij]. Having both 
heavy-ion and baseline p-\-p measurements with the same 
experiment is advantageous as it leads to a reduction of 
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the systematic uncertainties and, thus, to a more precise 
relative comparison of the data. 

In this paper, we present results on inclusive neutral 
pion production at mid-rapidity from proton-proton col- 
lisions at = 62.4 GeV from data collected during the 
RHIC 2006 run. A sample of events from longitudinally 
polarized p + p collisions (about 2 /3 of the total data 
sample) was used for double helicity asymmetry mea- 
surements. The other events from the 2006 data sam- 
ple were from transversely polarized p + p collisions and, 
along with the longitudinally polarized data, were used 
for the unpolarized cross section measurements, by aver- 
aging over the different initial spin states. 

The structure of this paper is as follows. The PHENIX 
subsystems used in this analysis are briefly introduced in 
section HIl The unpolarized tt'^ cross section analysis and 
the results are discussed in section Hill The tt'^ All anal- 
ysis and the results follow in section HVl and a summary 
is given in section IVl 



II. EXPERIMENT 

The PHENIX experiment at RHIC measured 7r°'s via 
7r° — > 77 decays using a highly segmented (A77 x Acj) ~ 
0.01 X 0.01) electromagnetic calorimeter (EMCal) [15| . 
covering a pseudorapidity range of \r]\ < 0.35 and az- 
imuthal angle range of A0 = tt. The EMCal comprises 
two calorimeter types: 6 sectors of lead scintillator sam- 
pling calorimeter (PbSc) and 2 sectors of a lead glass 
Cherenkov calorimeter (PbGl) . Each of the EMCal tow- 
ers was calibrated by the two-photon invariant mass from 
7r° decays and cross checked against the energy deposited 
by the minimum ionizing particles in the EMCal, and the 
correlation between the measured momenta of electron 
and positron tracks and the associated energy deposited 
in the EMCal. The uncertainty on the absolute energy 
scale was 1.2%. 

The 7r° data in this analysis were collected using two 
different triggers. One is a beam-beam counter (BBC) 
trigger which was defined by the coincidence of signals in 
two BBCs located at pseudorapidities ±(3.0 — 3.9) with 
full azimuthal coverage [16| . The time difference between 
the two BBCs was used to determine the collision vertex 
along the beam axis, which in this analysis was required 
to be within 30 cm from the center of the PHENIX inter- 
action region (IR). The other trigger is an EMCal-based 
high pt photon trigger, in which threshold discrimination 
corresponding to a deposited energy of 0.8 GeV was ap- 
plied independently to sums of analog signals from 2x2 
groupings of adjacent EMCal towers This trigger had 
limited efficiency for tt" detection at low pt (e.g. 50% 
in 1.0-1.5 GeV/c px bin) and close to 100% efficiency at 
PT >3 GeV/c. 

Beam-beam counters along with zero degree calorime- 
ters (ZDC) [13, which detect neutral particles near the 
beam axis {9 < 2.5 mrad), were utilized to determine 
the integrated luminosity for the analyzed data sample 



needed for the absolute normalization of the measured 
cross sections. Trigger counts defined with the BBCs and 
ZDCs were also used for the precise measurements of the 
relative luminosity between bunches with different spin 
configuration, and the spin dependence of very forward 
neutron production [1, [H, [3, detected by the ZDCs, 
served for monitoring the orientation of the beam polar- 
ization in the PHENIX interaction region (IR) through 
the run. These are necessary components of the spin 
asymmetry measurements. 



III. THE pp CROSS SECTION 

The unpolarized cross section analysis technique was 
very similar to our analyses of — 200 GeV data 0, 0] 
and is briefly discussed in section IIII Al Cross section 
measurements require an absolute determination of lu- 
minosity which is described in section IIIIBI The tt^ 
cross section results are presented and discussed in sec- 
tion ImCl 



A. n analysis 

The TT*^ yield in each px bin was determined from the 
two-photon invariant mass spectra. The background con- 
tribution under the tt*' peak in the two-photon invariant 
mass distribution varied from 75% in the lowest 0.5-0.75 
GeV/c Pt bin to less than 4% for pT > 3 GeV/c. 

One of the main corrections applied to the measured 
spectrum is the BBC trigger bias /^o , which is defined as 
the fraction of high pT tt^ events in the mid-rapidity spec- 
trometer acceptance which fire the BBC trigger. This 
fraction was determined from the ratio of the number of 
reconstructed tt" in the high px photon triggered sam- 
ple with and without the BBC trigger requirement. As 
shown in Fig.IU Jt^o was about 40% up to pt ~ 3 GeV/c 
and then monotonically dropped down to 25% a,t pT 1 
GeV/c. The drop can be explained by the fact that most 
of the energy is used for the production of high energy 
jets which contain the measured high pT tt*^ and there is 
not enough energy left to produce particles in the BBC 
acceptance 3.0 < I77I < 3.9, which was optimized for 
v/s — 200 GeV (where such a drop was not observed 
Q) and was not moved for the present ^/s = 62.4 GeV 
measurements. 

The main contributors to the systematic uncertain- 
ties of the measured tt^ spectrum are given in Table [H 
The "Energy scale" uncertainty includes uncertainties 
due to EMCal energy absolute calibration and nonlinear- 
ity. The "Yield extraction" uncertainty comes from the 
background subtraction. The "Yield correction" uncer- 
tainty comes from the correction for the geometric accep- 
tance, trigger efficiency, reconstruction efficiencies, detec- 
tor response, and photon conversion. The normalization 
uncertainty is not included and is discussed in section 

mm 
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FIG. 1: (color online) The fraction of the inclusive yield 
which satisfied the BBC trigger condition. 



TABLE I: Main systematic uncertainties in % of the tt" spec- 
trum from the PbSc for two representative pr bins (the PbGl 
uncertainties are similar). 



(pt) (GeV/c) 


1.2 


6.7 


Energy scale 


3.9 


13.1 


Yield extraction 


3.9 


2.0 


Yield correction 


6.4 


6.0 



The data sets from the two EMCal subsystems, PbSc 
and PbGl, were analyzed separately and combined for the 
final results. Results from the two subsystems were con- 
sistent within uncertainties. The systematic uncertainty 
of the combined result is reduced as the major systematic 
uncertainties in the two EMCal subsystems are not cor- 
related. For final 7r° cross section results, BBC-triggered 
events were used for px < 3 GeV/c and high pT pho- 
ton triggered events in coincidence with the BBC trigger 
were used for pr > 3 GeV/c. 



B. Vernier scan analysis 

The measured tt" cross section was normalized to the 
integrated luminosity for the analyzed data sample (L) 
which was determined from the number of BBC-triggered 
events using an absolute calibration of the BBC trigger 
cross section obbc- The value of gbbc is obtained via 
the van der Meer or Vernier scan technique [13]. This is 
a crucial part of the absolute cross section analysis and 
is therefore discussed in detail in this section. 

In a scan, the transverse widths of the beam overlap 
Ux and Uy were measured by sweeping one beam across 
the other in small steps while monitoring the BBC trigger 
rate. Then the instantaneous machine luminosity of each 



-'machine 



is computed as: 

In 



2TTax<7y 



■ Nb- 



Y, 



(3) 



where Nb and Ny are the bunch intensities of the two 
beams (^10^^/bunch), frev is the revolution frequency 
(78 kHz). The BBC trigger cross section ctbsc is the ra- 
tio of the BBC trigger rate when the beams were overlap- 
ping maximally (Rmax) to the effective luminosity L^ff'- 



O'BBC — Rmax/Leff, 



where 



L 



-'macriine 



^vertex; 



(4) 



(5) 



and Cvortex is the fraction of the number of collisions in the 
PHENIX interaction region (IR) within the BBC trigger 
vertex cut (usually \z\ < 30 cm). 

Lmachine was Corrected for the z dependence of the 
transverse beam sizes caused by the beam focusing in 
the IR (hour-glass effect) and for the beam crossing an- 
gle. The value of Cvortox was extracted from the ^-vertex 
distribution of events measured by the BBCs and was cor- 
rected for the dependence of the BBC trigger efhciency 
on the collision vertex position z along the beam axis. 
These corrections are discussed in more detail below. 

In p + p colhsions at ,/s = 62.4 GeV, the BBC trig- 
ger efficiency vs z shape was estimated from the com- 
parison with a "detector unbiased" ^-vertex distribution 
obtained from the convolution of colliding bunch inten- 
sity profiles along the z-axis as measured by Wall Cur- 
rent Monitors (WCMs) [2l[. The correction factor of 
0.83 ± 0.08 for Cvcrtox in Eq- ® was obtained, resulting 
in Cvcrtcx — 0.37 ± 10%. This approach is confirmed in 
p + p collisions at y/s — 200 GeV where the ZDCs have 
enough efficiency to measure the z vertex distribution. 
The efficiency of the ZDCs (located at z = ±18 m) does 
not depend on collision vertex position in the PHENIX 
IR, which was distributed with a sigma of 0.5 ~ 0.7 m 
around z = 0. The vertex distribution obtained with the 
WCMs is well reproduced by the measurement with the 
ZDCs at Vs = 200 GeV (Fig. Hi). 

Beam focusing in the IR causes bunch transverse sizes 
to vary away from the nominal collision point (z = 0) 
as a'^{z) = a^{z = 0) x (1 + z^/P*^), where P* is the 
value of the betatron amplitude function at the interac- 
tion point. This is the so-called hour-glass effect. The 
product (JxO'y in Eq. ^ should be replaced by an effec- 
tive {ax ■ cFy), which differs from what was measured in a 
scan (mainly due to the vertex cut implemented in BBC 
trigger). The correction due to this effect for Vernier 
scan data at ^/s = 62.4 GeV with a betatron amplitude 
function at the collision point of /3* = 3 m was simulated 
with WCM data and calculated to be 0.93 ± 0.02. The 
applicability of our calculational technique is illustrated 
in Fig. [2] with the high statistics Vernier scan data at 
= 200 GeV. 
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FIG. 2: (color online) Collision z-vertex distribution in the 
PHENIX IR measured by ZDCs in a Vernier scan at ^ = 200 
GeV (points) and calculations from convolution of colliding 
bunch intensity profiles along z-axis and including the hour- 
glass effect for /3* = 1 m, for bunches with typical length of 
1 m and transverse size of 0.3 mm (histograms); (a) beams 
are head-on; (b) one beam is 0.9 mm displaced relative to the 
other beam in the horizontal direction (illustrates the hour- 
glass effect) and (c) one beam is 0.9 mm displaced relative 
to the other beam in the vertical direction. The calculations 
include the bunch crossing angle with a vertical projection of 
0.15 mrad. 



pQCD ^=Pt 
MRST2002 PDF; fDSS FF 

NLL 
— NLO 




Figure [2b and [^t shows the sensitivity of our data for 
the transversely displaced beams to the hour-glass effect 
and to the crossing angle between the colliding beams, 
compared with a head-on vertex distribution in Fig. [2^. 
The two peaks in Fig. [2b and [5};, caused by the hour-glass 
effect, show an overlap of the diverging colliding beams 
at large \z\ in a particular displaced beam setting from a 
Vernier scan. The obvious asymmetry in the two peaks in 
Fig. [2t is a result of the non-zero crossing angle between 
colliding bunches. In all Vernier scan measurements the 
crossing angle was found to be less than 0.2 mrad, which 
translates to a negligible correction for Lmachine at ^/s = 
62.4 GeV, with a typical bunch length of 1 m and 
bunch transverse size of 1 mm. 

After all the corrections discussed above were applied, 
our BBC trigger cross section in p + p collisions at 
= 62.4 GeV was found to be (Jbbc = 13.7 mb with a 
systematic uncertainty of ±1.5 mb (±11%), i.e. ~ 40% 
of the world-average value of the inelastic p -\- p scat- 
tering cross section at -/s — 62.4 GeV[3|- Major con- 
tributors to the systematic uncertainty are 4% from the 
uncertainty in the normalization of bunch intensity mea- 
surements and in the calibration of the beam position 
measurements in the Vernier scan, 10% from the BBC 
trigger efficiency correction of evertex, and 2% from the 
hour-glass correction. 



C. TT cross section results and discussion 

Figure [3] presents the inclusive mid-rapidity tt^ invari- 
ant production cross section at -y/s = 62.4 GeV versus 
PT, from pT = 0.5 GeV/c to = 7 GeV [H. An over- 
all normalization uncertainty of 11% due to the uncer- 



FIG. 3; (color online) (a) The neutral pion production cross 
section at y/s = 62.4 GeV as a function of pt (circles) and 
the results of NLO (solid) and NLL (dashed) pQCD calcula- 
tions for the theory scale fi — pr- (b) The relative difference 
between the data and NLO pQCD calculations for the three 
theory scales jj, = pt/2 (upper line), pr (middle line) and 
2pT (lower line); experimental uncertainties (excluding the 
11% normalization uncertainty) are shown for the fj, = pr 
curve, (c) The same as b) but for NLL pQCD calculations. 



tainty in absolute normalization of the luminosity is not 
shown. The analyzed data sample with 0.76 x 10^ BBC 
triggers corresponded to about 55 nb^^ integrated lumi- 
nosity. The measurements fall within the large spread of 
ISR data [HEUl. 

The data are compared to NLO and NLL pQCD cal- 
culations at a theory scale fi = pt, where /i represents 
equal factorization, renormalization, and fragmentation 
scales [lH. The NLL corrections extend the NLO calcu- 
lations to include the resummation of extra "threshold" 
logarithmic terms which appear in the perturbative ex- 
pansion at not very high energies because the initial par- 
tons have just enough energy to produce the high pT par- 
ton that fragments into a final pion. The MRST2002 par- 
ton distribution functions [IJl and the fDSS set of frag- 
mentation functions (25j , which are extracted in NLO, 
are used in both NLO and NLL calculations. We have 
previously seen that the data are well described by NLO 
pQCD with a scale of /i = at Vs = 200 GeV [ii]. In 
contrast, NLO calculations with the same scale underesti- 
mate the tt" cross section at ^/s = 62.4 GeV. At the same 
time, it is known that NLO calculations are not always 
successful at describing low energy fixed target data 0, 
while NLL calculations have been successful [l3|- The 
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FIG. 4: (color online) (a) The neutral pion production cross 
section at ^/s = 62.4 GeV and -^s = 200 GeV as a func- 
tion of XT, scaled by (y^/GeV)" with n=6.38; the solid line 
is a parameterization of = 200 GeV data, (b) The pa- 
rameter n in ((6)l obtained from the ratio of invariant cross 
section at = 62.4 GeV and -^s = 200 GeV, at each 
XT of -^s — 62.4 GeV data; error bars show the statisti- 
cal and systematic uncertainties of the ^/s = 62.4 GeV and 
= 200 GeV data. The shaded band reflects the ll%e9.7% 
normalization uncertainty in the ^/s = 62.4 and 200 GeV 
cross section measurements, correspondingly. 



NLL calculations have a smaller scale dependence and 
describe our data well with /i = pT', however, as noted 
in [Ui subleading perturbative corrections to the NLL 
calculation may be significant. Therefore, the results 
may indicate that ^/s — 62.4 GeV is at an intermediate 
energy region where calculations that include threshold 
logarithm effects may describe the data more accurately. 
Therefore, below we show comparisons to both NLO and 
NLL at a scale of ^ = pt- 

General principles of hard scattering, including the 
principle of factorization of the reaction into parton dis- 
tribution functions for the protons, fragmentation func- 
tions for the scattered partons and a short-distance 
parton-parton hard scattering cross section, predicted a 
general xx-scaling form for the invariant cross section of 
inclusive particle production near mid-rapidity [2^ : 



1 



dp^ 



Pt 



-G{xt) 



(6) 



where xt = "^VtI \fs-, and F{xt) and G(xt) are universal 
functions. The parameter n relates to the form of the 
force-law between constituents. For example for QED 
or Vector Gluon exchange, n = 4 ^l^. Due to higher 
order effects, the running of the coupling constant q;((3^). 



the evolution of the parton distribution functions and 
fragmentation functions, and the initial-state transverse 
momentum /ct, n ys, not a constant but is a function of 
XT and n[xTT\fs) [11]. 

Figure 2^ shows the inclusive 7r° cross section scaled 
by \/s" for = 62.4 GeV and 200 GeV data 

[5| as a function of xt-, with the parameter n = 
6.38, which is a weighted average of n(xT) for xt > 
0.07 (corresponding to pt > 2 GeV/c at y/s — 
62.4 GeV). The parameter n{xT) was calculated as 
ln(fT62.4(a;T)/o'20o(a;T))/ ln(200/62.4) for eachxT of = 
62.4 GeV data, and ae2A and cr2oo are invariant differen- 
tial cross sections at y/i = 62.4 GeV and y^ = 200 GeV, 
respectively. Cross section values for the correspond- 
ing Xt at y^ = 200 GeV were obtained from pa- 
rameterization of the measured cross section at y/s = 
200 GeV: T(pt)(I+^4^ + " ^^P^))^' ^^P^) = 
1 where t = 4.5 GeV/c, w = 0.084 GeV/c, 

"3 po = L488 GeV/c, m = 10.82, 



l+exp((pT — ' 

A = 253.8 mb-GeV~2 

B = 14.7 mb-GeV-2+fe . c^-^, and k = 8.11. All 
y^ = 200 GeV data points agree with the parameter- 
ization curve within uncertainties. The parameterization 
is shown as the solid curve in Fig. 2^. 

At low XT, where soft physics dominates particle pro- 
duction, n{xT) is supposed to increase with xt due to 
the similar exponential shapes of the soft part of the in- 
variant cross section versus pt at different y^ (~ e~^^^) 
(27t . In the hard scattering region n{xT) is expected to 
decrease with increasing xt , due to stronger scale break- 
ing at lower pt- Such behavior of u^xt) is demonstrated 
by our data in Fig.[3)D. A similar drop in the parameter n 
at XT ^ 0.1 was observed at ISR energies [I4I. Figure |4|3 
also shows the possible transition from soft- to hard- 
scattering regions in ir'^ production at pt GeV/c. 
A similar conclusion was derived from the shape of the 
7r° spectrum at y^ = 200 GeV in This can serve as 
a basis for applying the pQCD formalism to the double 
helicity asymmetry data with pT > 2 GeV/c in order to 
allow access to AG. 



IV. INCLUSIVE tt" double HELICITY 
ASYMMETRY 

A. tt'' All analysis 

For the 2006 run, each of the two independent RHIC 
collider rings were filled with up to 111 bunches in a 
120 bunch pattern, with one of four, fiU-by-fiU alter- 
nating predetermined patterns of polarization sign for 
the bunches. Bunch polarization signs in each pattern 
were set in such a way that all four colliding bunch spin 
combinations occurred in sequences of four bunch cross- 
ings. That greatly reduced the systematic effects in spin 
asymmetry measurements due to variation of detector 
response versus time and due to possible correlation of 
detector performance with RHIC bunch structure. 
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To collect data from collisions of longitudinally polar- 
ized protons, the polarization orientation of the beams 
was rotated from vertical, the stable spin direction in 
RHIC, to longitudinal at the PHENIX IR and then back 
to vertical after the IR by spin rotators [29] . PHENIX 
local polarimeters measured the residual transverse com- 
ponent of the beam polarizations, using the spin depen- 
dence of very forward neutron production p, I TsI , 
observed by the ZDC, and by that means monitored the 
orientation of the beam polarization in the PHENIX IR 
throughout the run. 

The magnitudes of the beam polarizations at RHIC are 
measured using fast carbon target polarimeters (30] , nor- 
malized to absolute polarization measurements by a sepa- 
rate polarized atomic hydrogen jet polarimeter [31}. The 
luminosity- weighted beam polarizations over 11 RHIC 
fills used in the All analysis were (P)=0.48 for both 
beams, with 0.035 and 0.045 systematic uncertainty for 
the two RHIC beams, respectively. For the longitudi- 
nal polarization run period, the residual transverse po- 
larizations of the beams were (Pt/P)'^=0.11±0.15 and 
{Pt/PV 0.11 ± 0.12 for "Blue" and "Yellow" RHIC 
beams, respectively. The average transverse component 
of the product was (Pf • Pf )/(P^ • P^) < (Pt/P)^ ■ 
(Pt/P)^ = 0.012 ± 0.021; the average of the polariza- 
tion product over the run was (P^ • P^) — 0.23, with a 
systematic uncertainty of ±14%. 

Experimentally, the double helicity asymmetry for tt*^ 
production is determined as: 

LL |pB.pF| • + ^ 

where iV++ and iV_| are the number of tt'^'s and R is 

the relative luminosity between bunches with the same 
and opposite helicities. The analysis technique for the tt'^ 
All measurements is similar to our analyses of ^/s = 200 
GeV data [i,[li,[23. 

Double helicity asymmetry results were obtained from 
longitudinally polarized p + p collisions corresponding to 
~ 40 nb^^ integrated luminosity. Due to the limited 
BBC trigger efficiency for high px tt" events, high px 
photon triggered events without the BBC trigger condi- 
tion requirement were used for the asymmetry anal- 
ysis. This led to a slightly increased background in the 
7r° reconstruction and additional systematic uncertainty 
in the measurements of the relative luminosity between 
bunches with different helicity states. 

The background asymmetry under the tt^ peak in the 
two-photon mass distribution Af^^ was estimated from 
the counts outside the tt" peak, from a 177-217 MeV/c^ 
range in the two-photon mass distribution. Unlike our 
^/s = 200 GeV data analyses, a lower mass range was not 
used for Af^^ estimations due to cosmic background from 
non-collision events. This background contribution in the 
mass ranges of peak and higher mass was negligible 
(< 1%), and Af^^ was consistent with zero in all pt bins. 

Similar to our previous analyses, crossing-by-crossing 
accumulated number of BBC triggers were used for the 



measurements of the relative luminosity between bunches 
with different spin configuration. The uncertainty on the 
relative luminosity measurements SR was derived from 
the comparison between BBC trigger events and other 
trigger events, selecting different physics processes in dif- 
ferent kinematic ranges. In the y/s — 200 GeV data anal- 
ysis the comparison was done to triggers defined by the 
coincidence of signals from the two ZDCs 5, 18, 22]. Due 
to the limited efficiency of the ZDC at -^s = 62.4 GeV, 
the comparison in this analysis was performed with the 
number of events which fired simultaneously either of the 
two BBCs and either of the two ZDCs. Only 20% of the 
event statistics in this sample is contributed by BBC- 
triggered events, so this sample can be considered as es- 
sentially independent from the BBC event sample. From 
this comparison the upper limit of SR was estimated to 
be 0.6 X 10""^, which for the average beam polarizations of 
0.48 translates to SAll — 1.4 x 10""^, thepT independent 
uncertainty of the tt*^ double helicity asymmetry results. 
Single beam background < 0.35%, as determined by the 
trigger counts of non-colliding bunches and pileup prob- 
ability of < 0.02%, had negligible impact on the relative 
luminosity measurements. 

A transverse double spin asymmetry Axt, the trans- 
verse equivalent to Eq. ^ and ([7]), can contribute to All 
through the residual transverse component of the prod- 
uct of the beam polarizations discussed above. Similar 
to [5], Att was obtained from the sample with trans- 
verse polarization. The maximal possible Att effect on 
All was determined by ±5Att from the measured Att, 
which was < 0.15 • 6All in all pr bins. 



B. TT All results and discussion 

Figure [5] presents the measured double helicity asym- 
metry in tt'^ production versus pt [23| . A scale uncer- 
tainty of 14% in AJ^j^ due to the uncertainty in beam 
polarizations is not shown. The other systematic uncer- 
tainties are negligible, as discussed above, and checked 
using a technique to randomize the sign of bunch polar- 
ization and by varying the n*^ identification criteria [Tsj . 

Figure[5]also shows a set of curves from pQCD cal- 
culations that incorporates different scenarios for gluon 
polarization within the GRSV parameterization of the 
polarized parton distribution functions [s^. GRSV-std 
corresponds to the best fit to inclusive DIS data. The 
other three scenarios in FiglS] (GRSV-max, AG = 0, and 
AG = — G) are based on the best fit, but use the func- 
tions AG{xg) = G{xg),Q, —G{xg) at the initial scale for 
parton evolution {Q^ = 0.4 GeV^), where G{xg) is the 
unpolarized gluon distribution, and AG(a:g) is the dif- 
ference between the distribution of gluons with the same 
and opposite helicity to the parent proton. In Fig. O we 
compare our asymmetry data with both NLO and NLL 
calculations. The NLL calculations indicate that we have 
a reduced sensitivity to positive AG, but the effect is 
far less pronounced than at Fermilab fixed-target ener- 
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FIG. 5: (color online) The double helicity asymmetry for neu- 
tral pion production at -^i = 62.4 GeV as a function of pr 
(GeV/c). Error bars are statistical uncertainties, with the 
14% overall polarization uncertainty not shown; other exper- 
imental systematic uncertainties are negligible. Four GRSV 
theoretical calculations based on NLO pQCD (solid curves) 
and on NLL pQCD (dashed curves) are also shown for com- 
parison with the data (see text for details). Note that the 
AG = curves for NLO and NLL overlap. 



gies [ni- Similar to our = 200 GeV results [1, [11], 
our = 62.4 GeV All data do not support a large 
gluon polarization scenario, such as GRSV-max. 

Figure [S] presents the measured All versus xt in 7r° 
production overlaid with the results at y/s = 200 GeV Q . 
Clear statistical improvement can be seen at higher xt- 
For the measured range 2-4 GeV/c, the range of in 
each bin is broad and spans the range Xg = 0.06 — 0.4, as 
calculated by NLO pQCD . Thus our data set extends 
our Xg reach of sensitivity to AG and also overlaps pre- 
vious measurements, providing measurements with the 
same Xg but at a different scale. 



V. SUMMARY 

To summarize, we have presented the unpolarized cross 
section and double helicity asymmetries for tt" produc- 
tion at mid-rapidity, for proton-proton collisions at ^/s = 
62.4 GeV. The accuracy of the cross-section measure- 
ments, which fall within the large spread of ISR data, re- 
lies on direct 7r° two-photon decay reconstruction, precise 
calibration of the photon energy measurements, careful 
study of the trigger performance and accurate control of 
the integrated luminosity of the analyzed data sample. 
The results serve as a precise baseline for heavy-ion mea- 
surements. Comparisons to NLO and NLL theoretical 
calculations indicate that including the effects of thresh- 
old logarithms may be necessary to more accurately de- 
scribe the cross section at y/s = 62 A GeV. The All 
results extend the sensitivity to the polarized gluon dis- 
tribution in the proton to higher Xg compared to the 
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FIG. 6: (color online) The double helicity asymmetry for 
neutral pion production at = 62.4 GeV and 200 GeV 
as a function of xt- Error bars are statistical uncertain- 
ties, with the 14% (9.4%) overall polarization uncertainty for 
^ = 62.4 GeV (200 GeV) data not shown. Two GRSV the- 
oretical calculations based on NLO pQCD are also shown for 
comparison with the data (see text for details.) 



previous measurements at ^/s — 200 GeV. A prelimi- 
nary version of these double helicity asymmetry results 
was already used in a recent global fit of both RHIC and 
polarized DIS data to constrain AG [11]. 



Acknowledgments 

We thank the RHIC Polarimeter Group and the staff 
of the Collider-Accelerator and Physics Departments at 
Brookhaven National Laboratory and the staff of the 
other PHENIX participating institutions for their vital 
contributions. We thank W. Vogelsang and M. Strat- 
mann for providing the pQCD calculations and for in- 
formative discussions. We acknowledge support from the 
Office of Nuclear Physics in the Office of Science of the 
Department of Energy, the National Science Foundation, 
a sponsored research grant from Renaissance Technolo- 
gies LLC, Abilene Christian University Research Coun- 
cil, Research Foundation of SUNY, and Dean of the Col- 
lege of Arts and Sciences, Vanderbilt University (U.S. A), 
Ministry of Education, Culture, Sports, Science, and 
Technology and the Japan Society for the Promotion 
of Science (Japan), Conselho Nacional de Desenvolvi- 
mento Cienti'fico e Tecnologico and Fundagao de Am- 
paro a Pesquisa do Estado de Sao Paulo (Brazil), Nat- 
ural Science Foundation of China (People's Repubhc 
of China), Ministry of Education, Youth and Sports 
(Czech Republic), Centre National de la Recherche Sci- 
entifique. Commissariat a I'Energie Atomique, and In- 
stitut National de Physique Nucleaire et de Physique 
des Particules (France), Ministry of Industry, Science 
and Tekhnologies, Bundesministerium fur Bildung und 
Forschung, Deutscher Akademischer Austausch Dienst, 



10 



and Alexander von Humboldt Stiftung (Germany), Hun- 
garian National Science Fund, OTKA (Hungary), De- 
partment of Atomic Energy (India), Israel Science Foun- 
dation (Israel), Korea Research Foundation and Korea 
Science and Engineering Foundation (Korea), Ministry 
of Education and Science, Rassia Academy of Sciences, 



Federal Agency of Atomic Energy (Russia), VR and the 
Wallenberg Foundation (Sweden), the U.S. Civilian Re- 
search and Development Foundation for the Independent 
States of the Former Soviet Union, the US-Hungarian 
Fulbight Foundation for Educational Exchange, and the 
US-Israel Binational Science Foundation. 



[1] J. Ashman et al. (EMC), Phys. Lett. B206 364 (1988), [16 

Nucl. Phys. B328, 1 (1989). 
[2] A. Airapetian et al. (HERMES), Phys. Rev. D75, 012007 [17 

(2007); V.Yu. Alexakhin et al. (COMPASS), Phys. Lett. 

B647, 8 (2007). [18 
[3] G. Bunce et al, Ann. Rev. Nucl. Part. Sci. 50, 525 

(2000). [19 
[4] S.S. Adler et al. (PHENIX), Phys. Rev. Lett. 91, 241803 [20 

(2003) . 

[5] A. Adare et al. (PHENIX), Phys. Rev. D76, 051106 
(2007). 

[6] J. Adams et al. (STAR), Phys. Rev. Lett. 92, 171801 [21 

(2004) . 

[7] B.I. Abelev et al. (STAR), Phys. Rev. Lett. 97, 252001 [22 

(2006) . 

[8] S.S. Adler et al. (PHENIX), Phys. Rev. Lett. 98, 012002 [23 

(2007) . 

[9] P. Aurenche et al, Eur. Phys. J. C13, 347 (2000). [24 
[10] D. de Florian and W. Vogelsang, Phys. Rev. D71, 114004 [25 

(2005) . [26 
[11] D. de Florian, W. Vogelsang, and F. Wagner, Phys. 

Rev. D76, 094021 (2007). [27 
[12] A.L.S. Angelis (CCOR), Phys. Lett. B79, 505 (1978); 

C. Kourkoumelis et al. (R-806), Phys. Lett. B84, 271 [28 
(1979). 

[13] A.L.S. Angelis et al. (CMOR), Nucl. Phys. B327, 541 [29 
(1989); T. Akesson et al. (AFS), Sov. J. Nucl. Phys. 51, 
836 (1990); C. Kourkoumelis et al. (R-806), Z. Phys. C5, 
95 (1980); A. G. Clark et al. (CSZ), Phys. Lett. B74, 267 [30 
(1978); F. W. Busser et al. (CCRS), Nucl. Phys. B106, 
1 (1976); K. Eggert et al. (ACHM), Nucl. Phys. B98, [31 
49 (1975); F. W. Busser et al. (OCR), Phys. Lett. B46, [32 
471 (1973); F. W. Busser et al. (CCRS), Phys. Lett. B55, 
232 (1975); M. Banner et al. (Saclay), Nucl. Phys. B126, [33 
61 (1977); B. Alper et al. (BSC), Nucl. Phys. BlOO, 237 
(1975). 

[14] D. d'Enterria, J. Phys. G31, S491 (2005). [34 
[15] L. Aphecetche et al, Nucl. lustrum. Methods A499, 521 [35 
(2003). 



M. Allen et al, Nucl. lustrum. Methods A499, 549 

(2003) . 

C. Adler et al, Nucl. lustrum. Methods A470, 488 
(2001). 

S.S. Adler et al. (PHENIX), Phys. Rev. Lett. 93, 202002 

(2004) . 

Y. Fukao et al, Phys. Lett. B650, 325 (2007). 
S. van der Meer, ISR-PO/68-31, KEK68-64; A. Drees 
and Z. Xu, PAC-2001-RPAH116 Presented at IEEE Par- 
ticle Accelerator Conference (PAC2001), Chicago, Illi- 
nois, 18-22 Jun 2001. 

P. Cameron et al, Proc. of the Particle Accelerator Con- 
ference, p. 2146 (1999). 

S.S. Adler et al. (PHENIX), Phys. Rev. D73, 091102 
(2006). 

Tables of data available at http:/ /www.phenix.b nl.gov| 
/phenix /WWW /info /data/ppg087_dat a. html 

A. D. Martin et al, Eur. Phys. J. C 35, 325 (2004) 

D. de Florian et al, Phys. Rev. D75, 114010 (2007). 

R. Blankenbecler, S.J. Brodsky and J.F Gunion, Phys. 
Lett. B42, 461 (1972). 

S.M. Berman, J.D. Bjorken and J.B Kogut, Phys. 
Rev. D4, 3388 (1971). 

R. F. Cahalan, K. A. Geer, J. B. Kogut, and L. Susskind, 
Phys. Rev. Dll, 1199 (1975). 

W.W. MacKay et al, Proceedings of the 2003 Particle 
Accelerator Conference, edited by J. Chew, P. Lucas and 
S. Webber, IEEE, p.l697. 

O. Jinnouchi et al, RHIC/CAD Accelerator Physics 
Note 171 (2004). 

H. Okada et al, Phys. Lett. B638, 450 (2006). 

B. Jager et al, Phys. Rev. D67, 054005 (2003); M. Gliick 
et al, Phys. Rev. D63, 094005 (2001). 

M. Stratmann and W. Vogelsang, hep-ph/0702083 = J. 
Phys. Conf Ser. 69, 012035 (2007); W. Vogelsang, pri- 
vate communication. 

D. de Florian et al, Phys. Rev. Lett. 101, 072001 (2008). 
D.L. Adams et al, Phys. Lett. B 261, 197 (1991). 



